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ABSTRACT 

We present results from the high precision timing analysis of the pulsar-white dwarf 
(WD) binary PSR J1012-I-5307 using 15 years of multi-telescope data. Observations 
were performed regularly by the European Pulsar Timing Array (EFTA) network, 
consisting of Effelsberg, Jodrell Bank, Westerbork and Nangay. All the timing param- 
eters have been improved from the previously published values, most by an order of 
magnitude. In addition, a parallax measurement of tt = 1.2(3) mas is obtained for 
the first time for PSR J1012-I-5307, being consistent with the optical estimation from 
the WD companion. Combining improved 3D velocity information and models for the 
Galactic potential the complete evolutionary Galactic path of the system is obtained. 
A new intrinsic eccentricity upper limit of e < 8.4 x 10~^ is acquired, one of the 
smallest calculated for a binary system and a measurement of the variation of the 
projected semi-major axis also constrains the system's orbital orientation for the first 
time. It is shown that PSR J1012-)-5307 is an ideal laboratory for testing alternative 
theories of gravity. The measurement of the change of the orbital period of the system 
of Pb = 5(1) X 10""'^"' is used to set an upper limit on the dipole gravitational wave 
emission that is valid for a wide class of alternative theories of gravity. Moreover, it 
is shown that in combination with other binary pulsars PSR J1012-I-5307 is an ideal 
system to provide self-consistent, generic limits, based only on millisecond pulsar data, 
for the dipole radiation and the variation of the gravitational constant G. 
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1 INTRODUCTION 

PSR J1012-I-5307 is a 5.3 ms pulsar in a binary system 
with orbital period of 14.5 h and a low mass companion 
l|Nicastro et al.|[l995l ). It was discovered during a survey for 
short period pulsars with the 76 m Lovell radio telescope at 
Jodrell Bank. liTorimer et al.l (|l995l ) reported optical observa- 
tions revealing an optical counterpart within 0.2 ±0.5 arcsec 
of the pulsar timing position being consistent with a helium 
white dwarf (WD) companion. 

The optical observations of the WD companion provide 



* Member of the International Max Planck Research School (IM- 
PRS) for Astronomy and Astrophysics at the Universities of Bonn 
and Cologne 



unique information about the evolution of the binary system 
and the radio pulsar itself, such as comparing the cooling 
age of the companion with the spin-down age of the pulsa r 
(|Lorimer et al.lll995l : iDriebe et all 19981 : lErgma et aLll200ll ) . 

Using the N E2001 m odel for the Galactic distribution 
of free electrons (jCordes fc Lazio 2002) an d the pulsar's d is- 
persion measure (DM) of 9cm~'^ pc (|Nicastr"o et all 19951) a 
distan ce of ~ 410 pc is derived. In contrast, Callanan et al.l 
(1 19981 ) compared the measured optical luminosity of the WD 
to the value expected from WD models and calculated a dis- 
tance of d = 840 ± 90 pc. In addition they measured, by the 
Doppler shift of the measured H spectrum of the compan- 
ion, a radial velocity component of 44 ± 8kms~^ relative 
to the SSB. From the radial velocity and the orbital pa- 
rameters of the system the mass ratio of the pulsar and its 
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companion was measured to be g = mp/rric = 10.5 ± 0.5. 
Finally by fitting the spectrum of the WD to a grid of DA 
(hydrogen dominated) model atmospheres they derived a 
companion mass of rric = 0.16 ± 0.02 M©, a pulsar mass of 
rup = 1.64 ± 0.22 A'Iq and an orbital inclination angle of 
i = 52° ±4°. 

iLange et al.l (|200lh presented the most complete pre- 
cision timing analysis of PSR J1012+5307 using 4 years 
of timing data from the Effelsberg 100 m radio telescope 
and 7 years from the 76 m Lovell telescope. Using their low 
eccentricity binary model ELLl and combining the timing 
measurements with the results from the optical observations 
they derived the full 3D velocity information for the system. 
Furthermore, after correcting for Doppler effects, they de- 
rived the intrinsic spin parameters of the pulsar and a char- 
acteristic age of 8.6 ± 1.9 Gyr which is consistent with the 
WD age from the optical estimates. In addition, after calcu- 
lating upper limits for an extremely low orbital eccentricity 
they discussed evolutionary scenarios for the binary system 
but also presented tests and limits of alternative theories of 
gravitation. Finally, they discussed the prospects of future 
measurements of Post-Keplerian parameters (PK) which can 
contribute to the description of the orientation of the sys- 
tem and the calculation of stringent limits for the effective 
coupling strength of the scalar field to the pulsar. 

In this paper we revisit PSR J1012-I-5307 with seven 
more years of high-precision timing data and combined 
datasets from the European Pulsar Timing Array (EPTA) 
telescopes consisting of the 100 m Effelsberg radio-telescope 
of the Max-Planck-lnstitute for Radioastronomy, Germany, 
the 76 m Lovell radio-telescope at Jodrell Bank observa- 
tory of the University of Manchester, UK, the 94 m equiv- 
alent Westerbork Synthesis Radio Telescope (WSRT), the 
Netherlands and the 94 m equivalent Nangay decimetric Ra- 
dio Telescope (NRT), France. After a short description of 
the timing procedure and the technique of combining our 
multi-telescope data, we present the updated measurements 
of the astrometric, spin and binary parameters for the sys- 
tem. Specifically we show the improvement in all the timing 
parameters and in the orbital eccentricity limit and in ad- 
dition the value for the first time for PSR J1012-I-5307 of 
the timing parallax. Furthermore, we obtain a value for the 
orb ital period variati on, in agreement with the prediction 
of (|Lange et aLll200ll ). from which we test different theo- 
ries of gravitation and give one of the tightest bounds on 
a Parametrised Post-Newtonian (PPN) parameter. Finally, 
we present how the timing measurement of the change of 
the projected semi-major axis can complete the picture of 
the orientation of the binary system. 



2 OBSERVATIONS 
2.1 Effelsberg 

PSR J1012-I-5307 was observed regularly with the EflFels- 
berg 100 m radio telescope since October 1996 with typi- 
cal observing times of 5-15 min in three consecutive scans. 
Monthly observations were performed at 1400 MHz using 
the primary focus cooled HEMT receiver. It has a typi- 
cal system temperature of 25 K and an antenna gain of 
1.5KJy~^. In order to monitor dispersion measure (DM) 



variations, it was also observed irregularly until August 2006 
and monthly thereafter, at 2700 MHz. At these frequencies 
a cooled HEMT receiver located at the secondary focus was 
used which has a system temperature of 25 K. Finally, it was 
occasionally observed at 860 MHz using an uncooled HEMT 
receiver, located at the primary focus, with a typical system 
temperature of 60 K. The Effelsberg-Berkeley Pulsar Pro- 
cessor (EBPP) was used for coherent on-line de-dispersion 
of the signal from the LHC and RHC polarisations. It has 
32 channels for both polarisations spread across bandwidths 
of 40, 100 and 80 MHz a t 860, 1400 and 2700 MHz respec- 
tively (|Backer et al.|[l997l ). The output signals of each chan- 
nel were fed into de-disperser boards for coherent on-line de- 
dispersion and were synchronously folded with the topocen- 
tric period. 

Each TOA was obtained by cross-correlation of the pro- 
file with a synthetic template, which was constructed out of 
12 Gaussian com ponents fitted to a high sig nal-to-noise ratio 
standard profile jXramer et al.lll998l , ll999r ). The TOAs were 
locally time stamped using a H-maser clock at the observa- 
tory. They were converted to UTC using the GPS maser 
ofTset values measured at the observatory, and the GPS to 
UTC corrections were made from the Bureau International 
des Poids et Measures (BIPlvQ). 

2.2 Jodrell Bank 

PSR J1012-I-5307 has been observed with the Lovell radio 
telescope 2-3 times per month since 1993, at three diflerent 
frequencies. It is continuously observed at 1400 MHz and 
at 410 and 606 MHz, it was observed until 1997 and 1999, 
respectively. All the receivers are cryogenically cooled with 
system temperatures of 25, 50 and 35, respectively and their 
LHC and RHC polarisation signals are detected and inco- 
herently de-dispersed in a 2x32x0.0312 MHz filter bank at 
410 MHz, in a 2x6x0.1250 MHz filter-bank at 606 MHz and 
in a 2x32x1 MHz fiher-bank at 1400 MHz. The signals are 
synchronously folded at the topocentric pulsar period and 
finally copied to a disc. 

Each TOA was obtained by cross-correlation of the pro- 
file with a standard template, generated by the summation 
of high S/N profiles. The TOAs were transferred to GPS 
from a H-maser and the time stamp was derived as for Ef- 
felsberg. 

2.3 Westerbork 

PSR J1012-I-5307 was observed m onthly using the W SRT 
with the PuMa-I pulsar machine (jVoute et al.l 120021 '). We 
used three observing frequencies: observations at centre fre- 
quencies of 1380 MHz and 350 MHz were carried out monthly 
from August 1999, and the pulsar was observed occasionally 
at a centre frequency of 840 MHz from 2000 until 2002. The 
system temperatures were 27, 120 and 75 K, respectively and 
most observations were 30 minutes long. The WSRT obser- 
vations used a bandwidth of 8 x10 MHz for observations at 
840 MHz and 1380 MHz, and after September 2006 the 8 
bands were spread out over a total observing bandwidth of 
160 MHz for the 1380 MHz observations. The observations 
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at the low frequency setup used only two bands of 10 MHz, 
either centred at 328 and 382 MHz or 323 and 367 MHz. For 
the observations taken at 1380 or 840 MHz we used 64 fre- 
quency channels per 10 MHz band, and the observations at 
the low frequencies used 256 frequency channels per 10 MHz 
band. 

For each observation, the data were de-dispersed and 
folded offline. Integration over frequency and time resulted 
in one single profile for each observation. Each profile was 
cross-correlated with a standard template, generated by the 
summation of high S/N profiles, so finally only one TO A 
was computed for each observation. The TO As were trans- 
ferred to GPS from a H-maser clock and the time stamp was 
derived as for Effelsberg. 



2.4 Nan^ay 

PSR J1012+5307 was observed roughly every 3 to 4 weeks 
with the Nangay Radio Telescope (NRT) since late 2004. 
The NauQay Radio Telescope is equivalent to a 94-m dish, 
with a gain of 1.4KJy~^ and a minimal system tempera- 
ture of 35 K at 1.4 GHz in the direction of the pulsar. With 
the BON (Berkeley-Orleans-Nangay) coherent dedispersor, 
in the period covered by the observations, a 64 MHz band 
centred on 1398 MHz is split into sixteen 4 MHz channels 
and coherently dedispersed using a PC-cluster, with typical 
integration times of one hour. The Nangay data are recorded 
on a UTG(GPS) time scale marked at the analogue to digital 
converter by a Thunderbolt receiver (Trimble Inc.). Differ- 
ences between UTC and UTC(GPS) are less than 10 ns and 
therefore no laboratory clock corrections are needed. A sin- 
gle TOA was calculated from a cross-correlation with a pulse 
template for each observation of one hour. 



2.5 Multi-telescope precision timing 

Combing the EPTA multi-telescope datasets is not a triv- 
ial process. The main technique for achievin g the optimal 
combi nation of the data sets is presented bv ljanssen et al] 
l|2008h . In general, using different datasets from different 
telescopes and obtained at different frequencies requires ex- 
tra corrections, apart from the usual one of the transforma- 
tion of all the individual telescope arrival times to arrival 
times in the TAI at the solar system barycentre (SSB). The 
extra corrections needed are usually constant time offsets 
between different datasets of residuals. These offsets derive 
from differences in the procedure of calculating the TOAs at 
each telescope, specifically differences in the templates. The 
timmg software package TEMPC0 can fit for these time off- 
sets or "jumps" . In the current work seven of these jumps 
need to be fitted corresponding not only to the telescopes 
but also to the different frequencies used. Normally, three 
"jumps", one for each telescope, would be sufficient. How- 
ever, the TOAs at different frequencies are usually calcu- 
lated by different templates, which might not be aligned 
optimally. In the current case, this occurs for Effelsberg and 
Westerbork TOAs. In Table [T] the properties of the individ- 
ual datasets are presented. 
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The c ombination of the E PTA datasets has many ad- 
vantages (jjanssen et al.l l2008h . The need for continuous 
mult i- frequency TOAs for precise measurement of the dis- 
persion measure (DM) and monitoring of the DM variations 
was successfully accomplished. Most important the combi- 
nation of the high quality data from Effelsberg, Nangay, 
WSRT with the long time span data of Jodrell (& Effelsberg) 
provides us with a 15 year dataset of TOAs with no signifi- 
cant time gaps. Using all these EPTA datasets we improve 
and measure all the astrometric, spin and binary parameters 
of PSR J1012-I-5307 presented in the first column of Tabled 
As a comparison, in the second column the measured param- 
eters of only the current Effelsberg set is shown and in th e 
third the Effelsberg measurements from lLange et all (|200lh . 
From Table[5]it is clear that the EPTA provides the most ac- 
curate error estimations and in addition ~ 3(t measurements 
of 2 post-Keplerian parameters. 



3 ANALYSIS & RESULTS 

All the combined TOAs, weighted by their individual un- 
certainties, were analysed with TEMPO, using t he DE405 
ephemeris of the Jet Propu lsion Laboratory ( JPL) l|Standishl 
1 19981 . 120041 ') and the ELLl (|Lange et al.ll2001^ binary model. 
TEMPO minimises the sum of the weighted squared timing 
residuals, producing a set of improved pulsar parameters 
and the post-fit timing residuals. The uncertainties on the 
TOAs from each telescope are scaled by an appropriate fac- 
tor to achieve a uniform reduced — 1 for each data set. 
The best post-fit timing solution of all the combined resid- 
uals is presented in Figure [T] In the top panel the post-fit 
versus time is shown, with arbitrary offsets of the different 
datasets, where it is clear that the uncertainties of most of 
the data points are comparable. By comparing the parame- 
ters we get by different combinations of data sets (i.e. only 
the 1400 MHz data or without Jodrell) we concluded that it 
is much more efficient to finally use all the available datasets 
together, as shown in the lower panel of Figure [1] 



3.1 Timing parallax &: distance 

Apart from the common method of DM distance estimation 
(~ 410 pc) and the more rare optical one (d — 840 ± 90 pc) 
for PSR J1012-I-5307, there is another way of measuring the 
distance to a pulsar, with pulsar timing. In general, the tim- 
ing residuals of nearby pulsars demonstrate an annual par- 
allax. This timing parallax is obtained by measuring a time 
delay of the TOAs caused by the curvature of the emit- 
ted wavefronts at different positions of the Earth in its or- 
bit. This time delay has an amplitude of r|;g cos /3/(2cd) 



(|Lorimer fc Kramer |2005| ). where teq is the Earth-Sun dis- 
tance, /? the ecliptic latitude of the pulsar, c the speed 
of light and d the distance to the pulsar. This effect has 
been measured fo r very few pulsars like PSR B 1855-1-0 9 
iKaspi et all Il994l ), PSR J1713-H0747 llCamilo et al] Il994l ) , 
PSR J0437-4715 (ISandhu et all IT997I '). PS R J1744-1134 
dToscano et al.l Il999l ), PS R J2145-0750 (iLohmer et all 
l2004h and PSR J0030-H0451 (|Lommen et al.ll2006l ). Here for 
the first time we measure a parallax tt — 1.2 ± 0.3 mas for 
PSR J1012-I-5307. This parallax corresponds to a distance of 
d = 822 ± 178 pc which is consistent with the d = 840 ±90 pc 
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Table 1. Properties of the individual telescope datasets. 



Properties 


Effelsberg 


Jodrell Bank 


Westerbork 


Nangay 


Number of TOAs 


1972 


600 


234 


86 


Time span (MJD) 


50371-54717 


49221-54688 


51389-54638 


53309-54587 


Rms of individual set(fis) 


2.7 


8.6 


2.9 


1.9 


Observed frequencies (MHz) 


860, 1400, 2700 


410, 606, 1400 


330, 370, 840, 1380 


1400 



Table 2. Timing paranieters for PSR J1012+5307. Comparison between the datasets from EPTA, Effelsberg until 2008 and Effelsberg 
until 2001 l|Lange et al.l l2001^. 



Parameters 



EPTA 



Effelsberg 2008 



Effelsberg 2001 



Right ascension, a (J2000) 
Declination, 5 (J2000) 
l^a (masyr~i) 
fig (masyr^i) 
Parallax, w (mas) 

u (Hz) 
u (s-2) 
P (ms) 
P (ss-l) 

Epoch (MJD) 

Dispersion measure, DM (cm""^ pc) 

Orbital period, Pt, (days) 
Projected semi-major axis, x (It-s) 
r] (= esinLj) 
K (= e cos u>) 
Eccentricity, e * 

Longitude of the periastron, oj (deg) * 
Tasc (MJD) 

Pb (ss-i) 
X (ss"-') 

Solar system ephemeris model 

Number of TOAs 

RMS timing residual (/is) 



10'"12™33^4341010(99) 
53°07'02". 60070(13) 
2.562(14) 
-25.61(2) 
1.22(26) 

190.2678376220576(5) 
-6.20063(3)xl0-i« 
5.255749014115410(15) 
1.712794(9) X 10-20 

50700.0 

9.02314(7) 

0.60467271355(3) 
0.5818172(2) 
1.2(3)xl0-® 
0.06(31) X 10-6 
1.2(3)xl0-6 
93(14) 
50700.08162891(4) 

5.0(1.4) X 10-1'' 
2.3(8)xl0-i5 

DE405 
2892 
3.1 



10'"12™33''. 434089(13) 
53°07'02". 6001(2) 
2.56(2) 
-25.49(2) 
0.8(3) 

190.2678376220611(8) 
-6.20077(5) xlO-i'^ 
5.25574901411531(2) 
1.712833(13) x 10-20 

50700.0 

9.0209(3) 

0.60467271355(4) 
0.5818175(2) 
1.6(3)xl0-'^ 
0.14(34) xl0-'5 
1.6(3)xlO-'5 
85(12) 
50700.08162891(5) 

4(2)xl0-i"' 
< 1.8 X 10-1^ 

DE405 
1972 
2.7 



10'"12™33^ .43364(3) 
53°07'02". 5878(4) 
2.62(13) 
-25.0(2) 
< 1.6 

190.267837621910(3) 
-6.2070(5) X 10-" 
0.00525574901411947(7) 
1.71456(15)xl0-20 

50700.0 

9.022(3) 

0.6046727133(2) 
0.5818174(5) 
l.l(5)xl0-*^ 
0.20(50) X 10-6 
l.l(5)xl0-6 
79(24) 
50700.08162905(9) 

0.3(3)xl0-i2 
< 1.8 X 10-15 

DE200 
1213 
3.1 



* The eccentricity and the longitude of the periastron are calculated from the Laplace-Lagrange parameters, r] and k. 
Figures in parentheses are the nominal la TEMPO uncertainties in the least-significant digits quoted 



measured from the optical observations. The difference with 
the DM distance may point to a spars e free electron distri- 
bution in this location of the Galaxy (iGaensl er et al.l [20081 : 
IChatteriee et al. 2009) . By combining the optical and timing 
parall ax distance measurem ents we calculate the weighted 
mean (|Wall fc Jenkins! l2003l ) of the distance getting an im- 
proved value of d = 836 ± 80 pc. 



3.2 Improved 3D velocity measurement &c 
Galactic motion 

Combining the proper motion measured from timing (Table 
[2| and the distance to the system and the radial velocity 
of Vr = 44 ± 8 km-^s f rom the optical observations of the 
WD, iLange et al.l (|200ll ) managed to determine the full 3D 
motion of the pulsar relative to the SSB. Our new timing re- 



sults improve the proper motion measurements by an order 
of magnitude and using the combined parallax and optical 
distance we recalculate the 3D motion of the pulsar. We 
derive transverse velocities of 



Va ~ ^ad ~ 10.2 ± l.Okms^^ 
and 

^^,5 = fJ-sd — 101.5 ± 9.7kms~^ 



(1) 



(2) 



This yields a total transverse velocity of vt = 102.0 ± 
9.8kms~^. Using the radial velocity from the optical mea- 
surements we get the space velocity of the system Vspace ~ 
111.4 ± 9.5kms~^, consistent and almost three times more 
precise than the previous value. In addition, this value is 
still consistent with the average space velocity of millisec- 
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Figure 1. (Top) Post-fit timing residuals for the dataset of each 
telescope. From top to bottom, Effelsberg, Westerbork, Jodrell 
and Nangay. (Bottom) Post-fit timing residuals. Best timing so- 
lution with all the data sets yields the parameters in Table [2] 



Figure 2. (Top) Evolutionary path of PSR J1012+5307 on the 
Galactic plane. With arrow the current position of the pulsar and 
with * the position of the Sun is noted. (Bottom) The oscillations 
of PSR J1012+5307 above and below the Galactic plane through 
time. 



ond p ulsars of 130 km s ^ (|Lvne et al.lll99^ : iToscano et al.l 
1 19991 ). 

Since we know the 3D velocity of PSR J1012-I-5307 
we can try, for the first time, to track its Galactic path 
in time and space. Assuming a characteristic age of ~ 
10 Gyr and applying a model for the Ga lactic potential 
l|Kuiiken fc Gilmord Il989l : iPaczvriskil [l990l ) . we derive the 
evolutionary path of PSR J1012-f 5307 in the Galaxy from 
the point it started emitting as a millisecond pulsar. In the 
top part of Figure[21 the projection of the evolutionary path 
of the pulsar on the Galactic plane is shown, where the ar- 
row indicates the current position of the pulsar and the star 
indica tes the posit ion of the Sun, fo r the lKuiiken fc Gilmord 
l|l989h model (the lPaczvhskil l|l990l ) derives similar results). 
It is obvious that the pulsar is presently at one of its clos- 
est approaches to the Sun, which is why we can actually 
observe it. PSR J1012-f 5307reaches maximum distances of 
~ 30 kpc through its path, spending only a small fraction of 
its lifetime close to the solar system orbit. In Figure [5] (bot- 
tom) the movement of PSR J1012-I-5307 above and below 
the Galactic plane is shown versus time indicating that the 



pulsar is oscillating with a period of ~ 0.6 Gyr reaching a 
maximum distance of ~ 7 kpc above and below the Galactic 
plane. 

3.3 Eccentricity 

PSR J1012-f 5307 is a low eccentricity binary system. In our 
current timing solution we measure a value for t he eccentric- 
ity of (1.2 ± 0.3) X lO"*". However, as shown in lLange et al.l 
(|200lD . the Shapiro delay cannot be separated from the Roe- 
mer delay for this system, which leads to a small correction 
to this eccentricity value and specifically to the first Laplace- 
Lagrange parameter rj — esinoj. Thus, following their con- 
vention, for a companion mass of rric — 0.16(2) Mq, a mass 
ratio q — 10.5(5) and a mass function of fm ~ 0.000578 Mq 
we derive the range r and shape s of the Shapiro delay ac- 
cording to 

r = 4.9255 (mJMQ)^is (3) 
and 
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Figure 3. Distribution of values for the intrinsic eccentricity from 
Monte Carlo simulations. The dashed line cuts the distribution at 
the 95 per cent of the values. 



1/3 



(4) 



The intrinsic value of rj, calculated from equation (A22) of 
iLange et al.] HooJ), due to the contribution of the Shapiro 
delay, is = (— 1.4±3.4) x 10~^. The true eccentricity of the 
system is e = ^^/rf + K^ , where k — ecostj = (0.6 ±3.1) x 
lO"'^. By solving this equation in a Monte Carlo simulation 
(Figure[31), for datasets of the values and uncertainties of the 
intrinsic rj and k, we obtain an upper limit for the intrinsic 
eccentricity: 



e < 5.2 X 10" 
e < 8.4 X 10" 



(68 per cent C.L. 
(95 per cent C.L. 



(5) 
(6) 



This limit is bett er than the previously published value 
jLange et al.ll200ll ). 

This improved limit has another signifi cant impo r tance . 
Using the fluctuation-dissipation theorem, iPhinnevI (|l992l ) 
predicted that the orbital eccentricity, of a pulsar-WD bi- 
nary system is correlated with the orbital period. Specif- 
ically, there is the theoretical prediction of a relic orbital 
eccentricity due to convective eddy currents in the mass ac- 
cretion process of the neutro n star from the comp anion while 
in the red giant phase. In iLange et al.l (|200ll ) the eccen- 
tricity limit of PSR J1012-I-5307 was plotted versus the or- 
bital period and was compared with the model curves of the 
iPhinnev fc Kulkarnil ([l994) model. Our current eccent ricity 
limit is much lower than the one in iLange et all HooJ), but 
still in good agreement with the predictions from this model. 



3.4 Changes in the projected semi-major axis 

A change in the projected semi-major axis has been mea- 
sured in the current analysis, for the first time, for PSR 
1012+5307. The observed value of Xph s = 2.3(8) x 10"^'^ 
can b e the result of the various effects (|Lorimer fc Kramerl 
I2OO5I ): 

■D I -GW I deA -m , -so , -planet , -PM 
Xohs = X + X + + X + X + x^ +3; . (7) 



The first term, x^ , is the Doppler correction, which 
is th e combined effe ct of the proper motion of the sys- 
tem (Shklovskiil ll970l ) and a correction term for the Galac- 
tic acceleration. The contribution for the Galactic acceler- 



ation, X 



is of order 6 x 10 



Furthermore, we calcu- 



late the contribution of the Shklovskii effect to be x^^'' = 
x{fi^a + fJ-^)d/c ~ 8 x 10^^^. Both the contributions are very 
small compared to the observed value, thus, this term can 
be neglected. 

The second term, x'^^ is arising from the shrinking of 
the orbit due to gravitational-wave damping 

*°"--^f (I)" flip = 

(|Peterslll963 ). where T© = GMq/c^ = 4.9255 ^is and rUc is 
expressed in units of solar masses, x is the projected semi- 
major axis and Ph the orbital period. This contribution again 
is much smaller than the current measurement precision. 

The third term, deA/dt, is the contribution of the vary- 
ing aberration caused by geodetic precession of the pul- 
sar spi n axis, and is typically of order QP''°'^ P / Pf, « 2 x 
10"^* (|Damour fc Tavioij[l99a ). For a recycled pulsar, like 
PSR J1012-I-5307, the spin is expected to be close to parallel 
to the orbital angular momentum, which further suppresses 
this effect. Hence, the contribution is at least three orders 
of magnitude smaller than the observed. 

The fourth term, x^ , is representing a change in the size 
of the orbit caused by mass loss from the binary system. We 
investigate the mass loss due to the loss of rotation energy 
by the pulsar, which we consider as the dominant mass-loss 
effect. We initially calculate the change in the orbital period 
from the same contribution as follows: 



h = 8n 



Ir, 



4 p. 



10" 



(9) 



where M = + rric and the moment of inertia of the 
pulsar Ip ~ lO**^ gcm^. Subsequently, by Kepler's third law 
we calculate the change in the projected semi-major axis 
of the orbit to be ~ 10~^^. Thus, we can also neglect this 
contribution. 

The fifth and the sixth terms, x^° and i;P'<^"<=«, are the 
contributions due to the classical spin-orbit coupling caused 
by a spin-induced quadrupole moment of the companion and 
the existence of an additional planetary companion respec- 
tively. They can both be neglected. For the first one to be 
significant, a main -sequence star o r a rapidly rotatin g white 
dwarf companion (|Wex et al.lll998l : lKaspi et al.ll2000l) would 
be necessary. The second is not being considered because 
there is no evidence for another companion to the pulsar. 

Since all the other contributions are much smaller than 
the observed variation of the projected semi-major axis, we 
conclude that the measured value is arising from the last 
term of equation ([7|, i^^^. This is a variation of x caused 
by a change of the orbital inclinati on while the binary system 
is moving relat i vely to the SSB ([Arzoumanian et al.lll996l : 
iKopeikinI [l996l : ISandhu et al] 119971 ). The measurement of 
the effect is presented in the following equation: 



1.54 X 10 X cot i {—fid sin Q + fig cos Q) 



(10) 



where Q is the position angle of the ascending node. The 
quantities x, fia and ns are expressed in seconds and mil- 
liarcseconds per year, respectively. The proper motions and 
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Figure 4. Change of the projected semi-major axis versus po- 
sition angle of the ascending node. The two curves have been 
produced for i = 52° (peak at 180°) and i = 128° (peak at 0°). 
The sohd hne represents the measured value of x and the dashed 
lines the Irr and 2tr limits of x, from top to bottom. The latter 
constrains the orientation Q. 



the inclination angle have been measured and since we mea- 
sure the value of i^^^ = Xobs, we can, for the first time, 
restrict the orbital orientation Q of PSR J1012-I-5307. In 
Figure Uthe i^*^ versus the position angle of the ascending 
node is presented. Unfortunately, our measured value cannot 
fully restrict the orientation, however from the lower limits 
of Xobs we derive significant limits for the position angle. For 
an inclination angle of i = 52° we get 

151° <n< 220° (68 per cent C.L.) (11) 
and 

117° < n < 255° (95 per cent C.L. ), (12) 
while for i = 128° 

n<40° & n> 331° (68 per cent C.L.) (13) 
and 

n < 74° & n > 297° (95 per cent C.L. ). (14) 
3.5 Orbital period variations 

There are several effects that can contribute to changes in 
the observed orbital period of a binary system that can be 
either intrinsic to the orbit or just kinematic effects. The 
most important terms are: 



+ Pb 



(15) 



Pb, measured here for the first time 7\ = 5.0(1.4) x 10 
is the observable rate of change of the orbital period. The 
first and second terms, Pt , Pb , are the contributions 
from the mass loss from the binary and from tidal torques 
respectively. They can both be neglected in the case of 
PSR J1012-I-5307 because the first one is very small, as 
shown before (~ 10~^®), and the second is also small due 
to the lack of interaction between the pulsar and the com- 
panion. 



The third term, Pb , is identical to the first term of 
equation In order to account for the Galactic accelera- 
tion we have extended the lDamour fc Tavlorl l|l99lf ) expres- 
sion (for a flat rotation curve) to high Galactic latitudes 

Gal 



Kz \ sinbl 



cos I -f 



p 



+ sin^ I 



cos b , (16) 



where /3 = {d/ Rq) cosb — cos/. Kz is the vertical compo- 
nent o f Galactic acceleration taken from lHolmberg fc Flvnnl 
(2004), which for Galactic heights z = |dsin6| ^ 1.5 kpc can 
be approximated with sufficient accuracy by 



A-, (10" ■'cms"') ~2.27 2kpc + 3.68(l-e"*'^'"''P'^) , 

where zjjpc = 2;(kpc). Rq — 8.0=b0. 4 jEisenhauer et aT 
and^Q = 27.2±0.9kms-^kpc-i (|Feast fc Whitelock 



(17) 



2003) 



19971 ) 



are the Sun's Galactocentric distance and Galactic angular 
velocity (= Oort's A — B). For the pulsar's Galactic coordi- 
nates of I = 160.3° and b — 50.9° we find 



P^"'^ = (-5.6 ± 0.2) X 10" 



(18) 



We also calculate the contribution due to the Shklovskii 
effect according to the following: 



, Shk _ (nl + ^l)d 



Pb = (7.0 ± 0.7) X 10' 



(19) 



where we used the measured proper motion and the weighted 
mean of the distance discussed earlier, d. So, by summing 
we yield the Doppler correction: 

Shk 



+ Pb = (6.4 ±0.7) X 10" 



(20) 



The fourth term, Pb , is the contribution due to grav- 
itational wave emission. In general relativity, for circular or- 
bits it is given by 

■Gw -GR 192 7r /2^\5/3 (T^m^s^ 

^"^InJ (g-fi)V3 • (21) 



For PSR J1012-h5307 we find A = (-1.1 ± 0.2) x 10"". 

All the previous terms are the ones that are expected 
to contribute by using GR as our theory of gravity. How- 
ever, most alternative theories of gravity predict an extra 
contribution to the obse r ved o r bital period variation, via 
dipole radiation (see lWilll (jl993h : IWilll (|200ll ') and references 
therein). This dipolar gravitational radiation results from 
the difference in gravitational binding energy of the two bod- 
ies of a binary system, and is expected to be much larger 
than the quadrupolar contribution, especially if the binding 
energies of the two bodies of the binary system differ signif- 
icantly. Thus, the case of PSR J1012-I-5307 , where there is 
a pulsar- WD system, is ideal for testing the strength of such 
emission. One finds for small-eccentricity systems 

Pb = — 47r — — kdS , (22) 

Jo 

where rric is expressed in units of solar masses, ku refers to 
the dipole self-gravitational contribution, which takes dif- 
ferent values for different theories of gravity (zero for GR) 
and S = Sp — Sc i s the differe nce in the "sensitivities" of 
the two bodies (see I Willi (|l993l ) for definition), and /x is the 
reduced mass, mpiric/M, of the system. The sensitivity of a 
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body is related to its gravitational self-energy e. In the post- 
Newto nian limit s ~ e/mc^, which gives ^ 10"'* for a white 
dwarf (|Wili200ll ). Hence, we can neglect Sc in equation l|22|) 
since Sp ^ 0.2. Using the mass ratio q, equation (|22|l can be 
written as 



rydipole 



■ liDSp 



(23) 



, q + 1 

For a specific theory of gravity k_d is known and Sp can be 
calculated as a function of the equation-of-state of neutron 
star-matter. 

Finally, there are theories that predict that the locally 
measured gravitational constant G changes with time as the 
universe expands. A changing gravitational constant would 
cause a change in the orbital period, which for neutron star- 
WD systems can be written as 



G 
G 



(24) 



l|Damour et al.lll988l : lNordtvedtlll990h . 

The intrinsic change of the orbital period is the observed 
value minus the Doppler correction term from equation (|20|l : 

Pb"*" = Pb~ Pb = (-1.5 ± 1.5) X 10"" , (25) 
which agrees well with the GR prediction given above as 



(-0.4 ± 1.6) X 10" 



(26) 



Hence, there is no need for a P^^ or to explain the 
observed variation of the orbital period. On the other hand, 
this can be used to set limits for a wide class of alterna- 
tive theories of gravity, which we will show in the following 
sections. 



3.6 A generic limit for dipole radiation 

A tight system comprising a strongly self-gravitating neu- 
tron star and a weakly self-gravitating white dwarf should 
be a very efficient emitter of gravitational dipole radiation, 
if there is any deviation from general relativity that leads 
to a non- vanishing kd in equation p2|l . Hence, observations 
of such systems are ideal to constrain deviations of that 
kind. PSR J1012-I-5307 turns out to be a particularly useful 
system to conduct such a test, since: (1) the white-dwarf 
nature of the companion is affirmed optically, (2) the mass 
estimates in this double-line system are free of any explicit 
strong-field effect^, which are a priori unknown, if we do not 
want to restrict our analysis to specific theories of gravity, 
(3) the estimated mass of the pulsar seems to be rather high, 
which is important in the case of strong field effects that oc- 
cur only abo ve a certain critical mass, like the spontaneous 
scalarisation l|Damour fc Esposito-Fareselll993l ). 

In the previous section we have shown that the change 
in the orbital period is in full agreement with the prediction 
by general relativity, once the kinematic contributions are 
accounted for. Hence, any deviations from general relativity 
leading to a different P^^^ is either small or compensated for 



^ The mass estimation for the weakly self g ravitating white dwar f 
companion is done with Newtonian gravity llCallanan et al.ll998t) , 
and in any Lorentz-invariant theory of gravity the theoretical pre- 
diction for the mas s ratio does no t contain any explicit strong- 
field-gravity effects llDamoudbOOTi) . 



a potential 7^ . However, we can already limit the variation 
of the gravitational constant by using the published limit of 
G/G = (4±9) X 10~*^ yr~^ from the Lunar Laser Ranging 
(LLR) l| Williams et al.|[2004l ). In combination with equation 
pi)) it gives P^ = (— 1±3) x 10~*^, for the most conservative 
assumption Sp = 0. Hence, p^dipo'e ^ (_o,2 ± 1.6) x lO^^"*, 
which with the help of equation (|23p converts into 



KDSp = (0.5 ±6.0) x 10" 



(95 per cent C.L. 



(27) 



Furt hermore, if we assume Sp = Q.Hjnp/Mo) 
(c.f. iDamour fc Esposito-Faresd l|l992l )) we find 



K.D = (0.2 ± 2.4) X 10" 



(95 per cent C.L. 



(28) 



This number improves upon the previously published limit 
for PSR J1012-f5307 (Lange et al. 2001) by more than an 
order of magnitude. 

For the tens or-sca lar theories of 

iDamour fc Esposito-Fares3 l|l996f ) hldS^ ~ (op — ctcf' 
(op — ao)^ < 6 X 10~^, assuming that the efi^ective coupling 
strength of the companion WD to the scalar field, Oc, is 
much smaller than the pulsars and is approximately ao, 
where ao is a reference value of the coupling at infinity. This 
value im proves slightly on the previously published limit of 
7x lO"'^ (|Nice et al.l l2005'l. obtained from PSR J0751+1807. 
If the non-linear coupling parameter /3o is of order 10 or 
larger, then neutron stars are muc h more weakly couple d to 
the scalar field than white dwarfs (|Esposito-Faresell2005l ) . In 
this case, for PSR J1012-I-5307, (dp - Oc)^ « ao < 6 x 10"^ 
which is an order of magnitude wea ker than the limi t 
3.4 X 10"^ from PSR J1141-6545 l|Bhat et all 120081 '). 
Actually, in tensor-scalar theories of gravity the latter 
is possibly the most constraining pulsar binary system. 
However, since there has been no optical identification of 
the companion, that could establish its WD nature without 
mass determination based on a specific gravity theory, it 
is not yet possible to derive a general theory independent 
limit for dipole radiation from PSR J1141— 6545, as done 
here with PSR J1012-H5307. 

In the future, more accurate determination of the dis- 
tance and improvement of our Pt value, could further in- 
crease the precision of the PSR J1012-)-5307 limit. 



3.7 Combined limits on G and the dipole 
radiation with millisecond pulsars 

In the previous section we have used the LLR limit for G 
in order to provide a test for dipole radiation with a single 
binary pulsar system. On the other hand, a generic test for G 
cannot be done with a single binary pulsar, since in general 
theories that predict a variation of the gravitational constant 
typically also predict the existence of dipole radiation (|Willl 
ll993h PlFrom equations ((22)) and ((24| we can see that P^ oc 
Pi whereas pdipoio ^ p-i Y^ence, one can combine any two 
binary pulsars, with tight limits for Pi, and different orbital 
periods, in a joint analysis to break this degeneracy, and to 
provide a test for G and the dipole radiation that is based 
purely on pulsar data. A formally consistent way of doing 



■* It is interesting to point out, that in the Jordan-Fierz-Brans- 
Dicke theory p*?"'"^ _|_ pG _ g £qj. binary pulsars with white- 
dwarf companions that have orbital periods ~10 days. 
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this with white-dwarf binary pulsars is the application of 
equation 



ryexc 
f_b 

Pb 



1-1 



2M) 



(29) 



(see equations (I24|l and (|22|l ) to both binary pulsars, and 
solving in a Monte-Carlo simulation this set of two equa- 
tions for G/G and kd- This procedure properly accounts for 
the correlations due to this mutual dependence, and thus 
provides a self consistent test for G and the dipole radia- 
tion, that does not rely on LLR limits or theory specific 
assumptions. There remains the problem of getting a good 
estimate for Sp in a general theory independent test. As be- 
fore, we will use Sp = 0.1 [mp/Mo) keeping in mind that 
the limits given below are subject to certain changes, if a 
different assumption for Sp is made. 

With its short orbital period and its fairly well deter- 
mined masses PSR J1012-I-5307 is an ideal candidate for 
such a combined analysis. Presently, the bes t binary pulsar 
limit for G comes from PSR J0437-4715 (|Verbiest et all 
l2008l : iDeller et al.ll2008l ). where p^ipoio ^ q ^^^^ ^^^^ 
in the analysis to obtain the limit for G. Using this pulsar 
in combination with PSR J1012-f 5307 in a joint analysis as 
introduced above gives, with a 95 per cent C.L., 



G _ 
G " 
and 



(-0.7 ±3.3) X 10"^^ yr" 



(-0.009 ± 0.045) Jfo (30) 



(31) 



KD = (0.3 ±2.5) X 10"^ 

where Ho — 74kms~^ Mpc~^ has been used as a value 
for the Hubble constant (|Riess et al.ll2009l ). Our pulsar test 
therefore restricts G/G to less than a 20th of the expansion 
rate of the Universe. 

The li mit for G giv e n her e is clearly weaker than the 
one given in lDeller et al.l (|2008l ). The main reason for this is 
that the equation for used bv lDeller et al.l l|2008l ) does 
not account for the sensitivity of the pulsar as in equation 
l|24p. Furthermore, the combined analysis still allows for a 
certain range for p^ipoie pgj^ J0437-4715, leading to a 
somewhat weaker limit compared with an analysis that uses 
pdipoic _ seen in Fig. \S\ Although this limit 

for G is weaker than the LLR limit, it still provides a useful 
in dependent addi t ion to the LLR result, as has been argued 
in lVerbiest et al.l (120081 ). 

The limit for the dipole radiation is slightly weaker than 
the one given in the previous section. However, in contrast 
to the limit of the previous section, the limit here does not 
rely on the LLR result for G, and therefore constitutes an 
independent test based solely on binary pulsar observations. 

We would like to stress two facts about the advantage 
of combining specifically these two binary pulsars. Firstly, 
in both cases the companion white dwarf is identified opti- 
cally, and its non-compact nature is ascertained indepen- 
dently of the underlying theory of gravity. Secondly, the 
two pulsars seem to be rather heavy and similar in mass 
(~ 1.7M0 jfl, which is important in case we have efi'ects like 



^ In general, PSR J0437— 4715 does not allow the determi- 
nation of the pulsar mass, since this requires the mass func- 
tion, which contains explicit strong-field contributions. Within 



o 




2 4 

G/G [10-12 yr-i] 



Figure 5. Contour plots of the one and two a confidence regions 
on G/G and Kr> jointly. The elongation of the regions refiects 
the correlation due to the mutual dependence of the two systems, 
PSR J1012-I-5307 and PSR J0437-4715, in this combined test. 



spontaneous scalarisation above a critical neutron star mass 
(jPamour Sz EsTaosito-Fareselll993 ). In the future, more ac- 
curate measurements of Pt, and distance of the two pulsars 
could constrain even more our derived limits. 



4 CONCLUSIONS 

We have presented results from the high precision timing 
analysis of 15 yr of EPTA data for PSR J1012±5307. A 
first ever measurement of the timing parallax tt — 1.2(3) 
and distance has been obtained for this pulsar. Combined 
with information from optical observations of the WD com- 
panion an improved 3D velocity has been derived for the 
system. This information enables the derivation of the com- 
plete evolutionary path of the pulsar in the Galaxy, show- 
ing that it spent most of its lifetime far away from the so- 
lar system orbit. In addition, an improved limit on the ex- 
tremely low intrinsic eccentricity, e < 8.4 x 10~^ (95 per cent 
C.L.), has been acquired, which agrees well with the thcoret- 
ical e ccentricity-orbital period relation (Phi nney & Kulkarn| 
Il994h . 

Of particular interest is the measurement of the varia- 
tion of the projected semi-major axis, x — 2.3(8) x 10"^^ 
which is caused by a change in the orbital inclination as 
the system moves relative to the SSB. This measurement al- 
lowed us to set limits on the positional angle of the ascending 



the generic class o f con servative gravity theories l lWilll lll993l) : 
Id amour &: Taylor! lll992l )'). for instance, only the efli'ective grav- 
itational mass, Grrip, of PSR J0437— 4715 can be determined. 
However, if one assumes that Q deviates less than 20 per cent 
from G, the pulsar mass is in the range of 1.5 to 2.0 solar masses. 
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node, for the first time, tiie last unknown parameter in fully 
describing the orientation of this binary system. 

As a result of the significant measurement of the change 
in the orbital period of the system, 5.0(1.4) x 10"^"*, and 
the identified nature of the two bodies in this binary sys- 
tem, tests for alternative gravity theories could be per- 
formed. Firstly, a stringent, generic limit for the dipole ra- 
diation has been obtained from PSR J1012-I-5307, kds^ = 
0.5±6.0x 10"^ (95 per cent C.L.), with the use of the G limit 
from LLR. Secondly, in a self consistent analysis we have 
used PSR J1012-H5307 together with PSR J0437-4715 to 
derive a combined limit on the dipole radiation and the vari- 
ation of the gravitational constant, kd = (0.4 ± 2.6) x 10"'^ 
and G/G = (-0.7 ±3.3) x 10"^^ yr'^ (95 per cent C.L.) re- 
spectively. These limits have been derived just with the use 
of millisecond pulsar- WD binaries and are valid for a wide 
class of alternative theories of gravity. 
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